Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep.

BACKGROUND: The pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome. RESULTS: The 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained. CONCLUSION: According to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-ß/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

Expression, Polymorphism, and Potential Functional Sites of the BMPR1A Gene in the Sheep Horn.

Sheep horns are composed of bone and sheaths, and the BMPR1A gene is required for cartilage and osteogenic differentiation. Therefore, the BMPR1A gene may have a function related to the sheep horn, but its relationship with the sheep horn remains unclear. In this study, we first utilized RNA sequencing (RNA-seq) data to investigate the expression of the BMPR1A gene in different tissues and breeds of sheep. Second, whole-genome sequencing (WGS) data were used to explore the functional sites of the BMPR1A gene. Lastly, the allele-specific expression of the BMPR1A gene was explored. Our results indicate that BMPR1A gene expression is significantly higher in the normal horn groups than in the scurred groups. Importantly, this trend is consistent across several sheep breeds. Therefore, this finding suggests that the BMPR1A gene may be related to horn type. A total of 43 Single-Nucleotide Polymorphisms (SNPs) (F-statistics > 0.15) and 10 allele-specific expressions (ASEs) exhibited difference between the large and small horn populations. It is probable that these sites significantly impact the size of sheep horns. Compared to other polled species, we discovered ten amino acid sites that could influence horn presence. By combining RNA-seq and WGS functional loci results, we identified a functional site at position 40574836 on chromosome 25 that is both an SNP and exhibits allele-specific expression. In conclusion, we demonstrated that the BMPR1A gene is associated with horn type and identified some important functional sites which can be used as molecular markers in the breeding of sheep horns.

The SLC19A1-AS/miR-1343/WNT11 axis is a novel positive regulatory ceRNA network governing goat granulosa cell proliferation.

Long noncoding RNAs (lncRNAs), as competitive endogenous RNAs (ceRNAs), can directly or indirectly affect the proliferation and apoptosis of granulosa cells by regulating microRNA (miRNA) pathways. A ceRNA network of the SLC19A1-AS-miR-1343-WNT11 axis was constructed via comprehensive transcriptome sequencing of ovaries from goats with various fertility levels to further elucidate the function and regulatory mechanism of SLC19A1-AS in modulating miR-1343 and WNT11 during granulosa cell proliferation and apoptosis. Subsequent validation experiments were conducted in vitro using granulosa cells. In these experiments, we performed RNA immunoprecipitation (RIP) and identified SLC19A1-AS as a ceRNA in goat granulosa cells that promoted proliferation. Through bioinformatics prediction, luciferase reporter gene assays, and RNA pulldown assays, we confirmed that SLC19A1-AS acts as a sponge for miR-1343, preventing its binding to WNT11 mRNA and thereby increasing the expression of WNT11. This interaction also influenced the proliferation and apoptosis of granulosa cells. Our study systematically validated the biological function of the lncRNA-miRNA-mRNA ceRNA network in goat ovaries and revealed the potential regulatory mechanism by which SLC19A1-AS functions as a ceRNA in granulosa cells. These findings are expected to provide an important experimental foundation for further elucidating the physiological regulatory network of the ovary and contributing to reproductive health in goats.

Screening of Litter-Size-Associated SNPs in NOX4, PDE11A and GHR Genes of Sheep.

In previous studies, NOX4, PDE11A and GHR genes have been screened as important candidate genes for litter size in sheep by using the GWAS method; however, neither their effects on litter size nor the loci associated with litter size have been identified. In this study, three candidate loci (c.1057-4C > T in NOX4, c.1983C > T in PDE11A and c.1618C > T in GHR) were first screened based on our previous resequencing data of 10 sheep breeds. After the three loci were genotyped using Sequenom MassARRAY technology, we carried out population genetics analysis on the three loci and performed association analysis between the polymorphism of the three loci and the litter size of sheep. The results of population genetics analysis suggested that c.1057-4C > T in NOX4 and c.1983C > T in PDE11A may be subject to natural or artificial selection. The results of association analysis indicated that litter size was significantly associated with c.1057-4C > T in NOX4 and c.1983C > T in PDE11A (p < 0.05) in Small Tail Han sheep, and there was no significant interaction effect between the two loci on the litter size. In summary, c.1057-4C > T in NOX4 and c.1983 C > T in PDE11A can be considered candidate molecular markers for improving litter size in sheep.

A mutation in LPAR2 activates the miR-939-5p-LPAR2-PI3K/AKT axis to regulate the proliferation and apoptosis of granulosa cells in sheep.

Lysophosphatidic acid receptor-2 (LPAR2) is a G protein-coupled receptor, which is involved in various physiological processes such as cell development, proliferation, and apoptosis, and is thought to play an important role in follicular development and reproduction. There is evidence that miRNA recognition elements (MRE) in the gene 3'UTR often contain single nucleotide polymorphisms (SNPs) that can alter the binding affinity of the target miRNA, leading to dysregulation of gene expression. In this study, we detected a SNP in LPAR2 3 'UTR (rs410670692, c.*701C > T) in 384 small-tailed Han sheep using Sequenom MassARRAY®SNP genotyping. Association analysis showed that the SNP was significantly associated with litter size. Then, the effect of LPAR2 rs410670692 mutation on gene expression in sheep hosts was studied by molecular biotechnology. The results showed that the expression of LPAR2 in the TT genotype was significantly higher than that in the CC genotype, which confirmed the existence of rs410670692, a functional SNP, in LPAR2 3'UTR. We then used bioinformatics methods and double luciferase reporter gene assay to predict and confirm LPAR2 SNP rs410670692 as the direct targeting regulatory element of miR-939-5p. Cell transfection experiments further found that SNP rs410670692 down-regulated the mRNA and protein levels of LPAR2 by influencing the binding of miR-939-5p. To understand the function and mechanism of miR-939-5p in sheep granulosa cells (GCs), we conducted cell proliferation and apoptosis experiments which showed inhibited GCs proliferation along with promoted GCs apoptosis upon overexpression of miR-939-5p. Moreover, overexpression of miR-939-5p promotes apoptosis of granulosa cells by blocking the LPAR2-dependent PI3K/Akt signaling pathway. In conclusion, these results indicate that the SNP rs410670692 of LPAR2 is related to the litter size of small-tailed cold sheep, and miR-939-5p can act as a regulatory element binding to the C mutation of rs410670692 to regulate the expression of LPAR2, affect the development of GCs, and thus indirectly affect the litter size of sheep. These studies provide evidence for the involvement of LPAR2 polymorphism in sheep reproduction and are expected to provide new insights into the molecular genetic mechanisms of litter size traits in sheep.

Thyroid transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in small tail han sheep with FecB BB genotype.

BACKGROUND: The thyroid gland is an important endocrine gland in animals that secretes thyroid hormones and acts on various organs throughout the body. lncRNAs are long non-coding RNAs that play an important role in animal reproduction; however, there is a lack of understanding of their expression patterns and potential roles in the thyroid gland of the Small Tail Han (STH) sheep. In this study, we used RNA-Seq technology to examine the transcriptome expression pattern of the thyroid from the luteal phase (LP) and follicular phase (FP) of FecB BB (MM) STH sheep. RESULTS: We identified a total of 122 and 1287 differential expression lncRNAs (DELs) and differential expression mRNAs (DEGs), respectively, which were significantly differentially expressed. These DELs target genes and DEGs can be enriched in several signalling pathways related to the animal reproduction process. CONCLUSIONS: The expression profiles of DELs and DEGs in thyroid glands provide a more comprehensive resource for elucidating the reproductive regulatory mechanisms of STH sheep.

Detection of polymorphisms in six genes and their association analysis with litter size in sheep.

Litter size in sheep is a complex trait controlled by micro-effective polygenes. APAF1, CLSTN2, CTH, PLCB1, PLCB4, and CHST11 are all involved in mammalian reproduction. However, the effects of these genes on litter size in sheep are still unclear. Therefore, in this study, we used Sequenom MassARRAY® SNP assay technology to type the single nucleotide polymorphisms (SNPs) loci of six genes in five sheep breeds. The results showed that most sheep breeds contain three genotypes at each locus. Then, we conducted population genetic analysis on the SNPs of six genes and found that the polymorphic information content in all sheep breeds ranged from 0 to 0.37, and most sheep breeds were in Hardy-Weinberg equilibrium (p > 0.05). In addition, association analysis in Small Tail Han sheep indicated that the rs399534524 locus in CLSTN2 was highly associated with first parity litter size, and litter size in ewes with CT genotype was higher than that in ewes with CC genotype or TT genotype. Furthermore, the rs407142552 locus in CTH was highly associated with second parity litter size in Small Tail Han sheep, and litter size in ewes with CT genotype was higher than that in ewes with TT genotype. Finally, we predicted the CTH and CLSTN2 protein interaction network and found that HTR1E, NOM1, CCDC174 and ALPK3 interact with CLSTN2 and have been reported as candidate genes related to litter size in sheep. These results suggest that they may be useful genetic markers for increasing litter size in sheep.

Transcriptomic Analysis Reveals Differentially Expressed Circular RNAs Associated with Fecundity in the Sheep Hypothalamus with Different FecB Genotypes.

Circular RNAs (circRNAs) are a specific type of noncoding RNA, and some have defined roles in cellular and biological processes. However, little is known about the role of circRNAs in follicular development in sheep with FecB (fecundity Booroola) mutations. Here, the expression profiles of circRNAs were investigated using RNA sequencing (RNA-seq) in the follicular phase (F) and the luteal phase (L) of FecB mutant homozygous (BB) and wild-type (WW) Small Tail Han sheep. A total of 38,979 circRNAs were identified, and 314, 343, 336, and 296 of them were differentially expressed (DE) between BB_F and BB_L, WW_F and WW_L, BB_F and WW_F, and BB_L and WW_L, respectively. The length, type, and chromosome distribution of the circRNAs and the expression characteristic between the circRNAs and their host genes in the sheep hypothalamus were ascertained. Enrichment analysis showed that the host genes of DE circRNAs in the follicular and luteal phases were annotated to MAPK, gap junctions, progesterone-mediated oocyte maturation, oocyte meiosis, and other hormone-related signaling pathways, and the different FecB genotypes were annotated to the gap junctions, circadian entrainment, MAPK, and other hormone-related signaling pathways. The competing endogenous RNA network prediction revealed that the 129 target miRNAs might be bound to 336 DE circRNAs. oar_circ_0000523 and oar_circ_0028984, which were specifically expressed during the follicular phase in the BB genotype sheep, probably acted as miRNA sponges involved in the regulation of LH synthesis and secretion. This study reveals the expression profiles and characterization of circRNAs at two phases of follicular development considering different FecB genotypes, thereby providing an improved understanding of the roles of circRNAs in the sheep hypothalamus and their involvement in follicular development and ovulation.

Thyroid transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in Small Tail Han sheep with FecB++ genotype.

The thyroid gland is an important endocrine gland in animals, which mainly secretes thyroid hormones and acts on various organs of the body. Long-chain non-coding RNA (lncRNA) plays an important role in animal reproduction. However, there is still a lack of understanding of their expression patterns and potential roles in the thyroid of Small Tail Han (STH) sheep. In this study, RNA-seq was used to examine the transcriptome expression patterns of lncRNAs and mRNAs in the follicular phase (ww_FT) and luteal phase (ww_LT) in FecB++ genotype STH Sheep. A total of 17,217 lncRNAs and 39,112 mRNAs were identified including 96 differentially expressed lncRNAs (DELs) and 1054 differentially expressed mRNAs (DEGs). Functional analysis of genes with significant differences in expression level showed that these genes could be enriched in Ras signalling pathway, hedgehog (HH) signalling pathway, ATP-binding cassette (ABC) transporters and other signalling pathways related to animal reproduction. In addition, through correlation analysis for lncRNA-mRNA co-expression and network construction, we found that LNC_009115 and LNC_005796 trans target NIK-related kinase (NRK) and poly(A)-specific ribonuclease (PARN). LNC_007189 and LNC_002045 trans target progesterone-induced blocking factor 1 (PIBF1), LNC_009013 trans targets small mothers against decapentaplegic (SMAD1) are related to animal reproduction. These genes add new resources for elucidating the regulatory mechanisms of reproduction in sheep with different reproductive cycles of the FecB++ genotype STH sheep.

Expression profiles of oviductal mRNAs and lncRNAs in the follicular phase and luteal phase of sheep (Ovis aries) with 2 fecundity gene (FecB) genotypes.

FecB (also known as BMPR1B) is a crucial gene in sheep reproduction, which has a mutation (A746G) that was found to increase the ovulation rate and litter size. The FecB mutation is associated with reproductive endocrinology, such mutation can control external estrous characteristics and affect follicle-stimulating hormone during the estrous cycle. Previous researches showed that the FecB mutation can regulate the transcriptomic profiles in the reproductive-related tissues including hypothalamus, pituitary, and ovary during the estrous cycle of small-tailed Han (STH) sheep. However, little research has been reported on the correlation between FecB mutation and the estrous cycle in STH sheep oviduct. To investigate the coding and noncoding transcriptomic profiles involved in the estrous cycle and FecB in the sheep oviduct, RNA sequencing was performed to analyze the transcriptomic profiles of mRNAs and long noncoding RNAs (lncRNAs) in the oviduct during the estrous cycle of STH sheep with mutant (FecBBB) and wild-type (FecB++) genotypes. In total, 21,863 lncRNAs and 43,674 mRNAs were screened, the results showed that mRNAs had significantly higher expression levels than the lncRNAs, and the expression levels of these screened transcripts were lower in the follicular phase than they were in the luteal phase. Among them, the oviductal glycoprotein gene (OVGP1) had the highest expression level. In the comparison between the follicular and luteal phases, 57 differentially expressed (DE) lncRNAs and 637 DE mRNAs were detected, including FSTL5 mRNA and LNC_016628 lncRNA. In the comparison between the FecBBB and FecB++ genotypes, 26 DE lncRNAs and 421 DE mRNAs were detected, including EEF1D mRNA and LNC_006270 lncRNA. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functional enrichment analysis indicated that the DE mRNAs were enriched mainly in terms related to reproduction such as the tight junction, SAGA complex, ATP-binding cassette, nestin, and Hippo signaling pathway. The interaction network between DE lncRNAs and DE mRNAs indicated that LNC_018420 may be the key regulator in sheep oviduct. Together, our results can provide novel insights into the oviductal transcriptomic function against a FecB mutation background in sheep reproduction.

Association Analyses between Single Nucleotide Polymorphisms in ZFAT, FBN1, FAM184B Genes and Litter Size of Xinggao Mutton Sheep.

Previous studies have screened key candidate genes for litter size in sheep, including fibrillin-1 (FBN1), family with sequence similarity 184 member B (FAM184B) and zinc finger and AT-hook domain containing (ZFAT). Therefore, it is necessary to verify these genes in the Xinggao mutton sheep population and determine the associated loci for litter size. In this study, three loci (FBN1 g.160338382 T > C, FAM184B g.398531673 C > T and ZFAT g.20150315 C > T) were firstly screened based on the population differentiation coefficient between the polytocous and monotocous sheep groups. Then, population genetic analysis and association analysis were performed on these loci. The results revealed that the g.160338382 T > C in FBN1 was significantly associated with the litter size of sheep. Moreover, there was no significant interaction effect between the g.160338382 T > C locus and FecB on litter size. Notably, g.160338382 T > C is adjacent to the anterior border of exon 58 and belongs to a splice polypyrimidine tract variant, which may lead to alternative splicing and ultimately cause changes in the structure and function of the protein. In summary, our results provided a potentially effective genetic marker for improving the litter size of sheep.

Photoperiod Induces the Epigenetic Change of the GNAQ Gene in OVX+E2 Ewes.

GNAQ, a member of the alpha subunit encoding the q-like G protein, is a critical gene in cell signaling, and multiple studies have shown that upregulation of GNAQ gene expression ultimately inhibits the proliferation of gonadotropin-releasing hormone (GnRH) neurons and GnRH secretion, and ultimately affects mammalian reproduction. Photoperiod is a key inducer which plays an important role in gene expression regulation by affecting epigenetic modification. However, fewer studies have confirmed how photoperiod induces epigenetic modifications of the GNAQ gene. In this study, we examined the expression and epigenetic changes of GNAQ in the hypothalamus in ovariectomized and estradiol-treated (OVX+E2) sheep under three photoperiod treatments (short photoperiod treatment for 42 days, SP42; long photoperiod treatment for 42 days, LP42; 42 days of short photoperiod followed by 42 days of long photoperiod, SP-LP42). The results showed that the expression of GNAQ was significantly higher in SP-LP42 than in SP42 and LP42 (p < 0.05). Whole genome methylation sequencing (WGBS) results showed that there are multiple differentially methylated regions (DMRs) and loci between different groups of GNAQ. Among them, the DNA methylation level of DMRs at the CpG1 locus in SP42 was significantly higher than that of SP-LP42 (p < 0.01). Subsequently, we confirmed that the core promoter region of the GNAQ gene was located with 1100 to 1500 bp upstream, and the DNA methylation level of all eight CpG sites in SP42 was significantly higher than those in LP42 (p < 0.01), and significantly higher than those in SP-LP42 (p < 0.01), except site 2 and site 4 in the first sequencing fragment (p < 0.05) in the core promoter region. The expression of acetylated GNAQ histone H3 was significantly higher than that of the control group under three different photoperiods (p < 0.01); the acetylation level of sheep hypothalamic GNAQ genomic protein H3 was significantly lower under SP42 than under SP-LP42 (p < 0.05). This suggests that acetylated histone H3 binds to the core promoter region of the GNAQ gene, implying that GNAQ is epigenetically regulated by photoperiod through histone acetylation. In summary, the results suggest that photoperiod can induce DNA methylation in the core promoter region and histone acetylation in the promoter region of the GNAQ gene, and hypothesize that the two may be key factors in regulating the differential expression of GNAQ under different photoperiods, thus regulating the hypothalamus-pituitary-gonadal axis (HPGA) through the seasonal estrus in sheep. The results of this study will provide some new information to understand the function of epigenetic modifications in reproduction in sheep.

Thyroid Transcriptomics Revealed the Reproductive Regulation of miRNA in the Follicular and Luteal Phases in Small-Tail Han Sheep with Different FecB Genotypes.

MicroRNA (miRNA) is a type of endogenous short-stranded ncRNA that influences many biological processes such as animal growth, development and metabolism. The thyroid gland is an important endocrine gland in sheep, and an increasing number of studies have shown that the thyroid gland plays an important role in animal reproduction, but the molecular mechanisms of the thyroid gland in sheep reproduction are poorly understood. In this study, RNA-seq was used to detect transcriptome expression patterns in the thyroid gland between the follicular phase (FP) and luteal phase (LP) in FecB BB (MM) and FecB ++ (ww) small-tail Han (STH) sheep, respectively, and to identify differentially expressed miRNAs (DEMs) associated with reproduction. Bioinformatic analysis of the target genes of these DEMs revealed that they can be enriched in multiple GO terms associated with the reproductive process in animals and in the KEGG signaling pathway. The miRNA-mRNA coexpression network revealed that oar-miR-133 and oar-miR-370-3p may play an important role in sheep reproduction. The results of the dual-luciferase reporter assay suggest a possible targeting relationship between novel-51 and TARBP2. These results provided a novel resource for elucidating regulatory mechanisms underlying STH sheep prolificacy.

A body map of super-enhancers and their function in pig.

Introduction: Super-enhancers (SEs) are clusters of enhancers that act synergistically to drive the high-level expression of genes involved in cell identity and function. Although SEs have been extensively investigated in humans and mice, they have not been well characterized in pigs. Methods: Here, we identified 42,380 SEs in 14 pig tissues using chromatin immunoprecipitation sequencing, and statistics of its overall situation, studied the composition and characteristics of SE, and explored the influence of SEs characteristics on gene expression. Results: We observed that approximately 40% of normal enhancers (NEs) form SEs. Compared to NEs, we found that SEs were more likely to be enriched with an activated enhancer and show activated functions. Interestingly, SEs showed X chromosome depletion and short interspersed nuclear element enrichment, implying that SEs play an important role in sex traits and repeat evolution. Additionally, SE-associated genes exhibited higher expression levels and stronger conservation than NE-associated genes. However, genes with the largest SEs had higher expression levels than those with the smallest SEs, indicating that SE size may influence gene expression. Moreover, we observed a negative correlation between SE gene distance and gene expression, indicating that the proximity of SEs can affect gene activity. Gene ontology enrichment and motif analysis revealed that SEs have strong tissue-specific activity. For example, the CORO2B gene with a brain-specific SE shows strong brain-specific expression, and the phenylalanine hydroxylase gene with liver-specific SEs shows strong liver-specific expression. Discussion: In this study, we illustrated a body map of SEs and explored their functions in pigs, providing information on the composition and tissue-specific patterns of SEs. This study can serve as a valuable resource of gene regulatory and comparative analyses to the scientific community and provides a theoretical reference for genetic control mechanisms of important traits in pigs.

Genetic diversity, tissue-specific expression, and functional analysis of the ATP7A gene in sheep.

In humans, variation of the ATP7A gene may cause cranial exostosis, which is similar to "human horn," but the function of the ATP7A gene in sheep is still unknown. Tissue expression patterns and potential functional loci analysis of the ATP7A gene could help understand its function in sheep horn. In this study, we first identified tissue, sex, breed, and species-specific expression of the ATP7A gene in sheep based on the RNA-sequencing (RNA-seq) data. Second, the potential functional sites of the ATP7A gene were analyzed by using the whole genome sequencing (WGS) data of 99 sheep from 10 breeds. Last, the allele-specific expression of the ATP7A gene was explored. Our result showed the ATP7A gene has significantly higher expression in the big horn than in the small horn, and the ATP7A gene has high expression in the horn and skin, suggesting that this gene may be related to the horn. The PCA results show that the region around the ATP7A can distinguish horned and hornless groups to some extent, further indicating that the ATP7A may be related to horns. When compared with other species, we find seven ruminate specific amino acid sites of the ATP7A protein, which can be important to the ruminate horn. By analyzing WGS, we found 6 SNP sites with significant differences in frequency in horned and hornless populations, and most of these variants are present in the intron. But we still find some potential functional sites, including three missenses, three synonymous mutations, and four Indels. Finally, by combining the RNA-seq and WGS functional loci results, we find three mutations that showed allele-specific expression between big and small horns. This study shows that the ATP7A gene in sheep may be related to horn size, and several potential functional sites we identified here can be useful molecular markers for sheep horn breeding.

Comparative proteomics of ovaries elucidated the potential targets related to ovine prolificacy.

Small Tail Han (STH) sheep, a unique Chinese breed, is recognized for its early maturity, year-round estrus, and prolificacy. However, the molecular mechanism of its high prolificacy has not been fully elucidated. The Proteomics approach is feasible and effective to reveal the proteins involved in the complex physiological processes of any organism. Given this, we performed the protein expression profiling of ovarian tissues during the luteal phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PL vs. ML), and the follicular phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PF vs. MF), respectively. Parallel Reaction Monitoring (PRM) was conducted to validate the differentially abundant proteins (DAPs). The tandem mass tag (TMT) quantitative proteomic results showed that a total of 5,237 proteins were identified, of which 49 and 44 showed differential abundance in the PL vs. ML and PF vs. MF groups, respectively. Enrichments analyses indicated that the DAPs including TIA1 cytotoxic granule-associated RNA-binding protein-like 1 (TIAL1), nicotinamide phosphoribosyltransferase (NAMPT), and cellular retinoic acid-binding protein 1 (CRABP1) were enriched at the luteal phase, while TIAL1, inhibin beta-a-subunit (A2ICA4), and W5PG55 were enriched at the follicular phase, potentially mediating reproductive processes in polytocous ewes. Furthermore, six DAPs were verified using PRM, confirming the accuracy of the TMT data acquired in this study. Together, our work expanded the database of indigenous sheep breeds and provided new ovarian candidate molecular targets, which will help in the study of the genetic mechanisms of ovine prolificacy.

Comparative Transcriptomics Identify Key Pituitary Circular RNAs That Participate in Sheep (Ovis aries) Reproduction.

CircRNAs have been found to play key roles in many biological processes and have diverse biological functions. There have been studies on circRNAs in sheep pituitary, and some important circRNAs have been found. But there are still few studies on circRNAs in sheep pituitary with different fecundity. In this study, we obtained the circRNAs expression profiles in the pituitary of FecB ++ genotype Small Tail Han sheep with different fecundity and estrous phases. A total of 34,878 circRNAs were identified in 12 pituitary samples, 300 differentially expressed circRNAs (DE circRNAs) (down: 104; up: 196) were identified in polytocous sheep in the follicular phase (PF) and monotocous sheep in the follicular phase (MF) (PF vs. MF), and 347 DE circRNAs (down: 162; up: 185) were identified in polytocous sheep in the luteal phase (PL) and monotocous sheep in the luteal phase (ML) (PL vs. ML). Cortisol synthesis and secretion pathway (follicular phase) and estrogen signaling pathway (luteal phase) were obtained by functional enrichment analysis of circRNAs source genes. Competing endogenous RNA (ceRNA) network analysis of key DE circRNAs revealed that oar-circ-0022776 (source gene ITPR2, follicular phase) targeted oar-miR-432, oar-circ-0009003 (source gene ITPR1, luteal phase) and oar-circ-0003113 (source gene PLCB1, luteal phase) targeted oar-miR-370-3p. We also explored the coding ability of DE circRNAs. In conclusion, our study shows that changes in the pituitary circRNAs may be related to the response of the pituitary to steroid hormones and regulate the reproductive process of sheep by affecting the pituitary function. Results of this study provide some new information for understanding the functions of circRNAs and the fecundity of FecB ++ genotype sheep.

Mutation of the ETS1 3'UTR interacts with miR-216a-3p to regulate granulosa cell apoptosis in sheep.

ETS1, an important member of the ETS transcription factor family, is involved in a variety of physiological processes in living organisms, such as cell development, differentiation, proliferation and apoptosis, and is thought to be associated with embryonic development and reproduction. However, the polymorphism of ETS1 has been rarely studied, and its potential impact on the formation of reproductive traits in sheep remains unclear. Here, we first analyzed polymorphisms of ETS1 in a population of 382 small-tailed Han sheep with a lambing number record using the Kompetitive Allele Specific PCR (KASP) technique. The results showed the presence of a SNP locus rs161611767 (T > C) in the 3'UTR of ETS1. The association analysis showed the lambing number of first, second and third parity in the individuals with the CC genotype (2.51 ± 0.108, 2.51 ± 0.179, 1.27 ± 0.196) was higher than that of individuals with the TT genotype (1.79 ± 0.086, 1.56 ± 0.102, 0.56 ± 0.100) (P < 0.05). Then, molecular biotechnologies were used to investigate the effects of the EST1 rs161611767 mutant locus on host gene expression in sheep and the underlying mechanism of its effect on sheep reproduction. The RT‒qPCR results showed that the expression of ETS1 was higher in individuals with the CC genotype than in those with the TT genotype (P < 0.05). The dual luciferase reporter assay showed that the luciferase activity of ETS1 in sheep with the TT genotype was decreased compared to CC genotype (P < 0.05), confirming the existence of EST1 rs161611767 in the 3'UTR as a functional SNP. Given that the 3'UTR is an important regulatory region of gene transcription and translation, we performed bioinformatics prediction and confirmed that the SNP rs161611767 of ETS1 was a direct functional target of miR-216a-3p using dual luciferase activity assay, and the binding capacity of allele T was stronger than that of allele C. Subsequently, the cell transfection results showed that miR-216a-3p suppressed the endogenous expression of ETS1 in sheep primary granulosa cells (GCs). Finally, CCK-8, EdU, WB detection of marker proteins and flow cytometry were used to detect the effects of miR-216a-3p on GCs viability and proliferation/apoptosis, respectively. The results showed that miR-216a-3p inhibited the proliferation of GCs while promoting apoptosis of GCs. In conclusion, these results demonstrate that the SNP rs161611767 of ETS1 is associated with lambing number in small-tailed Han sheep, and miR-216a-3p can act as a regulatory element binding to the T mutation in rs161611767 to regulate ETS1 expression and affect GCs development, which may indirectly affect the number of lambs in sheep. These studies provide evidence for the involvement of ETS1 polymorphisms in sheep reproduction and are expected to provide new insights to elucidate the molecular genetic mechanisms of lambing traits in sheep.

Photoperiod Induces DNA Methylation Changes in the Melatonin Receptor 1A Gene in Ewes.

Research has shown that MTNR1A plays an essential role in the estrus cycle and seasonal reproduction changes in sheep. However, few people have focused on the DNA methylation of MTNR1A by season or photoperiod. In this study, using qPCR and Western blotting, we measured the MTNR1A expression in the hypothalamus of ovariectomized and estradiol-treated (OVX + E2) sheep under different photoperiod treatment conditions. Subsequently, the core promoter of the MTNR1A gene was identified, and its methylation level in sheep exposed to different photoperiod treatments was measured by pyrosequencing. The results showed that MTNR1A gene expression significantly differed between the short 42-day photoperiod (SP42) and the LP42 or combined SP-LP42 treatment groups (p < 0.05). In addition, we determined that the core MTNR1A promoter region was within 540 bp upstream of the transcriptional start site (TSS) and that the DNA methylation levels at CpG sites in the SP42 vs. LP42 and SP42 vs. SP-LP42 groups significantly differed. Finally, a significant negative correlation (p < 0.001) between gene expression and DNA methylation levels was revealed, suggesting that photoperiod may induce DNA methylation of the MTNR1A gene and thus change its expression. The findings provide valuable bases for the further study of seasonal reproduction in sheep.

Uterus proliferative period ceRNA network of Yunshang black goat reveals candidate genes on different kidding number trait.

Pregnancy loss that occurs in the uterus is an important and widespread problem in humans and farm animals and is also a key factor affecting the fecundity of livestock. Understanding the differences in the fecundity of goats may be helpful in guiding the breeding of goats with high fecundity. In this study, we performed RNA sequencing and bioinformatics analysis to study the uterus of Yunshang black goats with high and low fecundity in the proliferative period. We identified mRNAs, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) by analyzing the uterine transcriptomes. The target genes of the identified miRNAs and lncRNAs were predicted, and miRNA-mRNA interaction and competitive endogenous RNA (ceRNA) networks were constructed. By comparisons between low- and high-fecundity groups, we identified 1,674 differentially expressed mRNAs (914 were upregulated, and 760 were downregulated), 288 differentially expressed lncRNAs (149 were upregulated, and 139 were downregulated), and 17 differentially expressed miRNAs (4 were upregulated, and 13 were downregulated). In addition, 49 miRNA-mRNA pairs and 45 miRNA-lncRNA pairs were predicted in the interaction networks. We successfully constructed a ceRNA interaction network with 108 edges that contained 19 miRNAs, 11 mRNAs, and 73 lncRNAs. Five candidate genes (PLEKHA7, FAT2, FN1, SYK, and ITPR2) that were annotated as cell adhesion or calcium membrane channel protein were identified. Our results provide the overall expression profiles of mRNAs, lncRNAs, and miRNAs in the goat uterus during the proliferative period and are a valuable reference for studies into the mechanisms associated with the high fecundity, which may be helpful to guide goat to reduce pregnancy loss.